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PACS 87. 18. Gh - Cell-cell communication; collective behavior of motile cells 
PACS 64 . 60 . Cn - Order-disorder transformations 
PACS 64 . 70 . D- - Solid-liquid transitions 

Abstract - In this work we investigate the collective behavior of self-propelled particles that 
deform due to local pairwise interactions. We demonstrate that this deformation alone can induce 
alignment of the velocity vectors. The onset of collective motion is analyzed. Applying a Gaussian- 
core repulsion between the particles, we find a transition to disordered non-collective motion 
under compression. We here explain that this reflects the reentrant fluid behavior of the general 
Gaussian-core model now applied to a self-propelled system. Truncating the Gaussian potential 
can lead to cluster crystallization or more disordered cluster states. For intermediate values of the 
Gaussian-core potential we for the first time observe laning for deformable self-propelled particles. 
Finally, without the core potential, but including orientational noise, we connect our description 
to the Vicsek approach for self-propelled particles with nematic alignment interactions. 
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q Introduction. — The study of the onset of collective 
^notion in systems of self-propelled particles was pioneered 
t by Vicsek et al. 1 1 . In their model, point-like particles are 
^assumed to move with constant magnitude of velocity in 
^) two spatial dimensions. The angular orientation of the ve- 
^-H locity of each particle is subjected to stochastic noise. As 
H a crucial further ingredient, polar alignment rules for the 
velocity orientations of neighboring particles are explicitly 
OO superimposed. When decreasing the noise amplitude or in- 



spherical shape. Collective behavior within this model has 
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creasing the particle density, a macroscopic polar ordering 



J transition occurs for the velocity orientations. 
J* This transition corresponds to the onset of collective 
.^motion. It was studied in more detail later [2j|4], also for 
the case of nematic alignment between the velocity vectors 
H |5|6[ and for more complicated systems |7|8| . Furthermore, 
1 systems of self-propelled particles that feature finite non- 
zero volume were investigated [9|[l0] . In these studies it 
was demonstrated that steric interactions between rod-like 
particles can induce ordering. 

Here, we investigate the collective behavior in a system 
of self-propelled particles that can deform. For this pur- 
pose, we use an approach that was derived from symmetry 
arguments to model the behavior of a single self-propelled 
particle 11 12 . The description couples the velocity vec- 
tor of the particle to a deformation tensor. This tensor 
characterizes the symmetric deformations starting from a 



been studied under explicitly imposed global coupling 
or explicitly imposed pairwise alignment rules [14] . As an 
interesting feature, the breakdown of collective motion has 
been observed in the latter case under compression when 
an additional Gaussian repulsive potential was applied. 

In the following, we extend this model to the case in 
which pairwise interactions directly deform the particles. 
Imagine, for example, a scenario of soft self-propelled 
spheres that deform when squeezed together. Such ob- 
jects will be called soft deformable self-propelled particles. 
We present a mechanism of velocity alignment mediated 
through such deformations. This can result in macroscopic 
ordered collective motion. It should be stressed that ex- 
plicit alignment rules are not additionally imposed. 

We investigate the impact of this alignment mechanism 
on the onset of collective motion. Applying an anisotropic 
Gaussian-core potential between the particles, hexagonal 
translational order appears in the orientationally ordered 
collectively moving state. As was found for an isotropic 
Gaussian-core potential and an explicitly imposed align- 
ment mechanism 
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the state becomes disordered un- 
der compression. We now explain that this reflects the 
reentrant fluid behavior under compression generic for 
the Gaussian-core model. Truncating (and shifting) the 
Gaussian-core potential, we find variants of the disordered 
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state under compression, such as an orient at ionally disor- 
dered cluster crystal and a cluster fluid-like state. 

Polar orientational velocity order in the collectively 
moving state results from a sufficient strength of the 
Gaussian-core potential. We observe laning behavior at 
intermediate values of the potential strength. Without the 
Gaussian-core potential, point-like particles as in the Vic- 
sek model are considered. The implicit alignment mecha- 
nism then leads to nematic orientational order of the ve- 
locity vectors. Including orientational noise, we can link 
our description on soft deformable self-propelled particles 
to corresponding Vicsek models [5J6] . Breakdown of orien- 
tational order is observed under increasing noise strengths. 

The model. — Here, we consider soft deformable self- 
propelled particles. That is, pairwisely repulsive interac- 
tions will directly lead to deformations of the particles. As 
mentioned above, a simple example is given by soft self- 
propelled spherical particles that deform when squeezed 
together. Because of the assumed particle softness, a soft 
repulsive potential force will be applied. 

Our model uses the same variables as introduced before 
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a velocity vector and a second-rank deformation 
tensor Si for each particle i. Si is symmetric and traceless. 
It can be related to the lowest-order symmetric deviation 



for each particle from a spherical to an elliptic shape 15 



The particle number is N (i — 1...N). Using these 
variables together with the position of each particle r$, we 
consider the following coupled equations of motion in two 
spatial dimensions: 



Vi, 
7 V. 



v, 2 v, ; 



(1) 

aSi-Vi + ^fi, (2) 
k Si + b (1 - 2 Si : Si) (v iV< - -v?l) 

c(l-2S; rS^fif;- iffl). (3) 



Here, I is the unity matrix. We choose 7 > 0. Then, in 
eq. ([2]), the first two terms on the right enforce non-zero 
velocity values leading to self-propulsion in the stationary 
state. The first term on the right of eq. ([3| induces relax- 
ation back to the undeformed spherical state for n > 0. 
Direct coupling between velocities and deformations is in- 
cluded through the terms with the coefficients a and b. In 
contrast to refs. 
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we set a < and b > 0. Then, 
through these terms, the velocity tends to orient paral- 
lel to the long axis of the elliptic deformation for each 
particle. Single interactionless particles move ballistically 
on straight trajectories. Unless stated otherwise, we set 
7 = 0.5, k = 1, a = —1, and b = 0.5. 

The force acting on particle i mimics soft steric inter- 
actions between the particles. It is included in the velocity 
equation with the coefficient \i as a soft core-to-core repul- 
sive force. In contrast to ref. 14 , we include anisotropy 
in the steric interaction force due to elliptic deformations. 
The soft interaction potential U therefore contains the de- 
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Fig. 1: Implicit alignment interaction for soft deformable self- 
propelled particles. Two particles are injected from the bottom 
and collide. A time series of particle states (bottom to top) is 
shown for each panel (a)-(c). The smaller the value of 11 the 
more effective is the alignment interaction: (a) \i — 1, (b) 
Li = 0.01, and (c) \i — 0. States of particle separation smaller 
than 5cr are highlighted in brighter color. 



formation tensor, 



(4) 



where r^- = r 



Tj and a > is a characteristic interaction 
length, follows as = — Y^j^i dU/drij. The parameter 
K is set to 1 here since its effect is already included by 
H (and c), and we choose a = 1. In realistic systems of 
self-propelled particles, the repulsive interaction potential 
is typically not long-ranged. We will therefore truncate 
the potential at a cut-off distance r c (and shift it). 

In the previous studies, the term with the coefficient c 
in eq. ^ has not been considered. It is the key part of this 
work. We introduce this term to model soft deformable 
self-propelled particles. Then the particles directly deform 
due to soft repulsive steric interparticle interactions. For 
this purpose, we set c > 0. At each time step, the defor- 
mation tensor for each particle can be parameterized as 
Si = Si(hifii — 1 1). Si denotes the degree of deformation, 
the unit vector hi marks the long axis of deformation. 
Si = corresponds to the undeformed spherical ground 
state and Si = 2 to a completely squeezed state. Here, we 
confine ourselves to medium deformations and limit the 
degree of deformation Si to the interval Si G [0, 1] through 
the prefactor (1 — 2S^ : S^). We stress again that there 
is no explicit alignment rule between the particle orienta- 
tions as it was considered, e.g., in refs. 13 14 . 



Orientational alignment mechanism. — First we 
explain the implicit alignment mechanism for soft self- 
propelled particles illustrated in fig. [l] The results were 
obtained by numerically integrating eqs. ([l])-([3| forward 
in time for c = 20. In the figure, deformation is indicated 
by the ellipsoidal shape of the particles. For this pur- 
pose, we parameterize the undeformed spherical bound- 
ary of particle i by Ro^(^) = i2o(cos(<p), sin(y?)) with Rq 
the radius and ip G [0,27r[. The boundary location of the 
deformed state is drawn via R^(^) = (I + ckS^) • Ro;(</?) 
with an extra factor of a = 2 used for illustration only 
in the figures. Orientation and magnitude of the veloc- 
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Fig. 2: Snapshots from the ordering process of a representa- 
tive sample, starting from a disordered state of random initial 
conditions. The sample is still in the disordered state in (a) 
at t = 600, but has orientationally and translationally ordered 
in (b) at t = 12500. System parameters: N = 300, p = 0.14, 
M = 0.01, c = 2.5. 
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Fig. 3: Time t Vop> o.q until a polar orientational degree of 
velocity order vop of 0.9 is reached as a function of the de- 
formational force strength c. Each data point is an average 
over 100 samples starting from a different disordered state of 
random initial conditions. The inset shows a typical time evo- 
lution of the sample averaged vop (here for c = 0.8). System 
parameters: N 140, p 0.24, p 0.05. 



ity vectors are marked by the orientation and length of 
the short thick lines within the ellipsoids, respectively. To 
draw the figures, we adjust the ground state radius Rq of 
the undeformed particles for best visualization. 

When two particles come close to each other ( "collide" ) , 
they deform due to steric interactions modeled by the term 
with the coefficient c in eq. ([3|. For c > the particles 
are compressed along the connecting line between their 
centers of mass. Consequently the particles elongate per- 
pendicularly to that direction. Since a < and b > 0, 
the velocity vector and long axis of deformation tend to 
align parallel to each other for each particle. Together, 
these two mechanisms lead to an effective alignment of 
the velocity vectors of the two particles. The smaller the 
repulsive force in eq. ([2|, the more efficient is the align- 
ment mechanism as illustrated in figs. [I] (a)-(c). 

Polar orientational alignment. — To investigate 
the collective behavior of soft deformable self-propelled 
particles, we first study the ordering process in an ini- 
tially disordered crowd. For this purpose, N particles were 
simulated in a periodic box of size L x L. In the numer- 
ical calculations, to obtain box lengths of order unity, all 
lengths were rescaled by a factor of 100<j. A measure for 
the density is defined by p = N<j 2 /L 2 . As for all subse- 
quently reported results, we started from a random spatial 
distribution, random values for the velocity orientations 
and magnitudes, and vanishing deformation. Snapshots 
from the ordering process for a representative sample can 
be found in fig. [2] for the given material parameters (only 
part of the numerical calculation box is shown). 

It is important to note that the alignment mechanism in 
eq. (|3| through the force term with the coefficient c alone 
would lead to nematic orientational alignment. From fig. [2] 
(b), however, we find polar orientational alignment of the 
velocity vectors. This polar alignment must therefore re- 
sult from the core repulsion in eq. ([2| with coefficient p. 
Also the hexagonal translational order must be attributed 
to the core repulsion. Both conclusions will be verified. 



To analyze the role of the implicit alignment mechanism, 
the time it takes until the system orientationally orders 
was measured as a function of c. For this purpose, N = 
140 particles were simulated in a periodic box with density 
p = 0.24. Other system parameters were the same as used 
to obtain fig. [I] except for p = 0.05. 

As an order parameter, we evaluated the magni- 
tude of the sample- averaged velocity orientation, vop = 
II ~k Y^iLi v i/^i||- We consider the sample as being ordered 
when vop > 0.9. Fig. [3] illustrates as a function of c the 
time r Vo P >o.9 that it takes for the sample to orientation- 
ally order. Each data point is an average over 100 samples 
of different initial conditions. It is obvious that alignment 
through deformation strongly supports the onset of collec- 
tive motion. The ordering time r VOP> o.g seems to diverge 
for small values of c for the given parameter values. In 
the inset of the figure we show the ordering process as a 
function of time for one specific value of c. 

The ordering process can also be tracked by the time 
evolution of the mean square displacement. An example 
is given by fig. [4] in a log- log plot. Initially, before the first 
collisions occur, we find a slope close to 2 indicating bal- 
listic motion. There follows a time interval of slope closer 
to 1 when particle velocities are not yet orientationally 
ordered. Therefore collisions are frequent, which leads to 
an overall behavior close to diffusion. We observe a kink 
in the curve when collective motion sets in. Finally, all 
particles move into the same direction and collisions are 
rare. Consequently the motion of each particle is close to 
ballistic again. As a result, the slope of the mean square 
displacement asymptotically approaches a value of 2. 

In a previous study on deformable self-propelled parti- 
cles, it was found that under compression the hexagonal 
translational and polar orientational order breaks down 
when a critical density is reached 
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The order reap- 
peared under subsequent expansion. A hysteretic effect 
was numerically observed. As an interparticle potential 
the isotropic analog of eq. Q was used. 

We investigated this effect in detail and found that de- 
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Fig. 4: Time evolution of the mean square displacement (msd) 
plotted on double logarithmic scales reflects the ordering pro- 
cess. The evolution started from random initial conditions. In 
the disordered state the slope of the msd data curve decreases 
(a line of slope 1 is included for comparison) . Long-ranged or- 
dering starts around the point where the curvature becomes 
positive. Asymptotically, the slope approaches 2 in the or- 
dered collectively moving state. The inset shows the degree of 
polar orientational velocity order vop for comparison. System 
parameters as given in the caption of fig. [2] 



formability of the particles is not necessary to observe it. 
If we fix the particle shape, we can still observe the break- 
down of order under compression. This is the case even for 
spherical particles. However, increasing the deformational 
interactions via higher values of c supports the effect, and 
the breakdown of order occurs at lower densities. 

As a consequence, we identify the soft interaction poten- 
tial of Gaussian form as the main source of the reported ef- 
fect. In the equilibrium case, Gaussian-core systems have 
been studied in detail [I6HI8], also for anisotropic forms 



of the Gaussian interaction potential 19 and anisotropic 
particles [20]. At high density, the corresponding phase 
diagram features melting of crystalline phases into a fluid. 
Since the system is typically fluid at low densities as well, 
this high-density fluidity is called reentrant. We therefore 
conclude that we observe the analog of this reentrant fluid 
behavior in self-propelled Gaussian-core systems. 

In fig. [5] we qualitatively show examples of the behav- 
ior that we found at high densities. We varied the cut-off 
distances r c of the Gaussian-core potential. At the initial 
density p = 0.2, all samples showed hexagonal polar or- 
der analogous to fig. [2](b). Example trajectories of single 
particles are indicated in white. The local line thickness 
decreases with increasing magnitude of the recorded veloc- 
ity. Only part of the numerical calculation box is shown. 

Figs. [5] (a) and (b) correspond to a cut-off distance 
r c = 5<j. After global order has broken down, we still find 
some local order (a), which vanishes under further com- 
pression (b). The velocity magnitudes along a trajectory 
in (b) are highly non- homogeneous, indicating frequent 




Fig. 5: Illustration of the collective behavior at high densi- 
ties after sufficiently slow compression from a hexagonally or- 
dered polar state. The initial density was p = 0.2 in all cases. 
Trajectories of selected particles are indicated in white (line 
thickness decreases with increasing velocity). Different pan- 
els correspond to different cut-off radii r c and densities: (a) 
r c = 5a, p = 0.84, (b) r c = 5a, p = 0.93, (c) r c = 2a, p = 0.76, 
(d) r c = 2a, p = 1.27, (e) r c = 1.5a, p = 1.27, (f) r c = 2.5a, 
p 0.37. System parameters: N 1000, p 0.05, c 2.5. 



collisions. This scenario is identified with the reentrant 
fluid behavior outlined above for the Gaussian-core model. 

Interesting behavior occurs at smaller cut-off lengths. 
Figs. [5] (c) and (d) show an example for r c = 2a. Again 
the softness of the potential plays an important role. It 
allows putting several particles at the same spatial loca- 
tion with only finite energetic penalty. Under slow com- 
pression, this happens already when the system is still in 
the ordered crystalline state (c) . Several lattice points are 
occupied by more than one particle. The motion is still 
orientationally ordered as indicated by the straight tra- 
jectories. At higher densities this order breaks down. We 

here in a system of self- 



obtain a cluster crystal 21 22 



propulsion, with several particles stacked at each lattice 
site (d). The deformed particles are not ordered orien- 
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Fig. 6: Laning in a system of soft deformable self-propelled 
particles of high aspect ratio. Red particles move downward, 
turquoise ones upward. Only part of the numerical calculation 
box is shown. System parameters: N = 1200, p = 0.14, a = 
-10, p 0.03, c= 1.2. 



tationally. As indicated by the localized trajectories, the 
particles are trapped at the lattice points and the clusters 
are immobile. The situation may be connected to the ef- 
fect that particles in a plane under compression stack on 
top of each other at high densities. Due to this effect, the 
area density p can become larger than 1. 

At still smaller cut-off lengths r c = 1.5<r, this structure 
is more mobile, as indicated by the trajectories in fig. [5] 
(e). In effect, the particles can move over the plane. We 
may call this scenario a "cluster fluid-like" state. 

A special situation was found for r c = 2.5<r. Before 
melting, the hexagonal crystal structure turned into a rect- 
angular one. The latter is still orientationally ordered, as 
can be seen from the straight trajectories in fig. [5](f). 

Slowly expanding the systems to their initial density, all 
of these structures turned back into the polar hexagonal 
crystal structure. We also observed hysteretic effects [14] . 

Laning. — Laning has not been observed previously 
for deformable self-propelled particles. For rigid self- 
propelled particles, it has been reported very recently 
23 24] . The laning state is characterized by globally ne- 
matic, but locally polar alignment of the velocity vectors. 
We have already seen that high values of the strength p of 
the Gaussian-core repulsion lead to global polar order. In 
the next section, we will see that without the Gaussian- 
core repulsion, p = 0, nematic order is obtained globally 
and locally. Thus the laning state should be found for 
intermediate values of the core repulsion strength p. In 
addition, the system parameters need to be chosen to give 
a high aspect ratio of the typical particle deformation. 

At this stage, we only point out that laning can indeed 
be observed under these circumstances. Fig. [6] shows an 
example. Again we started from random initial conditions 
to obtain this result. The velocity alignment is nematic on 
a global scale and polar within each lane. Nonzero but not 
too high values of c > seem to support the macroscopic 
separation into different lanes. In the example of fig.[6j we 
measured a sample-averaged degree of deformation (s) « 
0.8. This corresponds to an aspect ratio of about 2.33 : 1. 




Fig. 7: Nematic and polar orientational velocity order parame- 
ters vop and orientational deformational order parameter sop 
as a function of orientational deformational noise strength A r . 
Left and right insets illustrate nematic alignment below and 
orientational disorder above the order-disorder transition, re- 
spectively. Red particles move downward, turquoise ones up- 
ward; only a very small portion of the sample is depicted. Sys- 
tem parameters: N — 800, p — 0.09, a — —5, p — 0, c — 0.5. 



Nematic orientational alignment. — Finally, our 
scope is to connect our description of soft deformable self- 
propelled particles with implicit alignment mechanism to 
the Vicsek model. As a first step, we set the steric core-to- 
core repulsion to zero (/i = 0) . We thus consider sterically 
point-like particles. The only remaining alignment mech- 
anism via the term with the coefficient c is of nematic 
character. It should lead to nematic orientational order of 
the particle velocities. Indeed we observe such a scenario 
as shown, for example, by the left inset in fig. [7] 

Second, we have to introduce orientational noise into 
our equations. It turns out that noise in the orientation 
of the deformation tensors is much more effective than 
noise in the orientation of the velocity vectors v^. There- 
fore, at each time step and for each particle i we pick a 
random angle & from a Gaussian distribution and multi- 
ply it by a constant factor A r . The deformation tensor 
is then rotated by the angle A r &. In this way we introduce 
purely orientational noise, without changing the degree of 
deformation. A r sets the noise strength. 

An example plot of the orientational order parameters 
as a function of the noise strength A r is shown in fig. [7J As 
expected from the Vicsek model with nematic alignment 
interactions [5j[6], the nematic orientational velocity order 
parameter goes to zero at a critical noise strength. At 
least for the finite system sizes investigated, the transition 
appears to be continuous. However, it is not possible to 
draw a final conclusion about the nature of the transition 
because a discontinuity may emerge at very high particle 
numbers 3,6 . At all noise strengths, the degree of defor- 
mational order is practically of the same magnitude as the 
degree of nematic velocity order. As expected, the polar 
velocity order is close to zero for all values of A r . We have 
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thus linked the scenario of soft deformable self-propelled 
particles to the corresponding Vicsek approach. 

Conclusions. — We introduced a description of soft 
deformable self-propelled particles. For this purpose, we 
included into the equations of motion a pairwise repul- 
sive force that can directly deform the particle shape. As 
demonstrated this leads to an implicit pairwise velocity 
alignment mechanism. The latter can significantly accel- 
erate the process of long-ranged orientational ordering. 

Furthermore, we studied the breakdown of collective 
motion under compression for the case of the introduced 
anisotropic repulsive Gaussian-core potential. This tran- 
sition is identified as the analog to the reentrant fluid be- 
havior in Gaussian-core equilibrium systems [l6 19 . The 
new soft deformational interactions were observed to shift 
this transition to lower densities. Interesting scenarios 
such as rectangular translational order, cluster crystals, 
and less ordered cluster states were observed at small cut- 
off lengths for the Gaussian interaction potential. We are 
not aware of any previous studies on cluster crystals or 
crystals of rectangular lattice structure for self-propelled 
particles. In addition, we for the first time found laning 
in deformable self-propelled particle systems. 

When strong Gaussian-core repulsion enforces transla- 
tional crystalline order, only polar orientational velocity 
order is compatible with the imposed lattice structure. At 
very low values of the core repulsion, nematic velocity or- 
der emerged from our implicit alignment mechanism. Lan- 
ing was found for intermediate values as a state of global 
nematic but local polar orientational velocity order. 

Finally, we linked our characterization to the Vicsek 
approach for point-like particles with nematic alignment 
interactions. Increasing the noise of deformation orien- 
tations caused the nematic orientational order to break 
down at a critical noise strength. 

Our approach should be useful for the investigation of 
soft biological tissue. As an example we mention soft 
motile cells that migrate on flat substrates [25]. Possi- 
ble additional adhesive forces can be included into our 
description. Another example may be given by soft vesi- 
cles that deform but do not merge when approaching each 
other. Candidates for non-biological model systems are 
self-propelled droplets at gas-liquid interfaces 26 or in 
bulk fluid 27 . These droplets were demonstrated to de- 
form under steric confinement interactions [26] and to sur- 
vive collisions 27 . Finally, we wish to point out the qual- 
itative similarity in appearance between our alignment 
mechanism illustrated in fig. [I] and the process of inelas- 
tic collisions in general. The latter was recently outlined 
as a bridge connecting particle-based descriptions of the 
Vicsek- type to the characterization of granular media [28] . 
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